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Abstract. QED radiative corrections have been calculated for leptonic and hadronic variables in parity- 
violating elastic ep scattering. For the first time, the calculation of the asymmetry in the elastic radiative 
tail is performed without the peaking-approximation assumption in hadronic variables configuration. A 
comparison with the PV-A4 data validates our approach. This method has been also used to evaluate 
the radiative corrections to the parity-violating asymmetry measured in the GO experiment. The results 
obtained are here presented. 

PACS. 25.30.Bf Elastic electron scattering ~ 13.60.-r Photon and charged-lepton interactions with hadrons 
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1 Introduction 

Elastic scattering of longitudinally polarized electrons is 
subject to parity violation through the interference be- 
tween 7 and Z° exchange. These experiments give ac- 
cess to the weak nucleon form factors (FF) , which are the 
equivalent, in the weak sector, of the usual electromagnetic 
form factors Ge and Gm- The weak nucleon form factors 
are related in turn to the strange form factors G^ and 

G , which are the contributions of strange currents to 
the form factors (see jT| and the following review arti- 
cles [21EE1EI)- According to QCD, this strangeness con- 
tribution arises from the presence of ss pairs in the nucleon 
sea. Many experiments have been performed recently or 
are stiU running at Bates (SAMPLE ^2,7,8 ), Mainz (PV- 
A4 miinj) and Jefferson Lab (G° [Hj and HAPPEX 12 ). 

Electromagnetic radiation produced from the emission 
of a real or virtual photon by the electron (incoming or 
outgoing) or by the target (before or after interaction), 
gives rises to a radiative tail which extends to very low 
energies (in theory, down to zero energy for the scattered 
electron). Since detectors have an experimental resolution 
and since cuts are used in the data analysis, the mea- 
sured cross section and asymmetry have to be corrected 
in order to be compared to theoretical models. The first 
calculations applied to elastic ep scattering were done by 
Tsai ^21 ) followed by a series of review papers [nill5lll6| . 
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This formalism has been later extended to scattering of 
polarized electrons 17 . All these calculations were done 
for experiments in which the scattered electron is detected, 
which was the case of SAMPLE, PV-A4, HAPPEX or G° 
at backward angles ^Hj ■ The originality of the G^ experi- 
ment at forward angles was the detection of recoil protons. 
In this case, some of the approximations commonly used 
when scattered electrons are detected, such as the peak- 
ing approximation, are no more valid. Thanks to its large 
mass, radiative emission from the proton is negligible but 
the proton kinematics is affected by the radiative emission 
from the electron (angle, energy, Q^). 



QED radiative corrections have been calculated for 
hadronic kinematic variables in ep elastic scattering |19| 
and applied to recoil proton polarization. In this case, a 
method based on an electron structure-function represen- 
tation, which is the analog of the Drell-Yann representa- 
tion [201, was used. These calculations were applied to ep 
scattering experiments done at Jefferson Lab |21| . aim- 
ing to determine the ratio of electric to magnetic proton 
form factors G^/G^at high momentum transfer as pro- 
posed by Akhiezer and Rekalo The classical method 
for computing corrections is based on the separation of 
the momentum phase space into hard- and soft-photons 
contributions to avoid infrared divergences This in- 
troduces a cutoff parameter which makes this method not 
easily applicable to construct an event operator. Bardin 
and Shumeiko 23 proposed a covariant cancellation me- 
thod of infrared divergences which does not introduce ad- 
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ditional parameters, however a cutoff has still to be intro- 
duced for generating real radiated photons. 

In the present work we calculate the corrections in 
both the leptonic and hadronic variables using an origi- 
nal method which is free of infrared divergences. Its in- 
terest is that it is exact and it can be integrated easily 
into numerical simulation programs such as GEANT. We 
then apply it to G° forward angles using the GO-GEANT 
code Pl]. Another original feature of the present calcu- 
lation is the inclusion of exchange, in addition to 7 
exchange, allowing to calculate the electroweak asymme- 
try of the radiative tail. We then calculate the corrected 
quantities of interest: integrated number of counts, time- 
of-flight (tof) spectra, distributions. A full account of 
the present calculations is given in the thesis of Hayko 
Guler gS] (in French). 

In section |5] we develop the formalism and define the 
Lagrangians, in section|3we describe our method for avoid- 
ing divergences, in section^lwe describe the apparatus 
and simulation method. The results are given in section 
and we conclude in sectional 



2 Theoretical formalism for elastic scattering 

The data analysis of parity- violating electron-nucleon scat- 
tering experiments involve the extraction of an asymmetry 
in the helicity-correlated cross section. The raw data are 
first converted into an experimentally measured asymme- 
try (Aeip). That means that the false asymmetry due to 
helicity-correlated fluctuations in intensity, energy, posi- 
tions and angles of the electron beam have already been 
taken into account. We assume also that background sub- 
traction has already been done. 

The asymmetry is commonly defined as: 



A ^ 



I da -y I da 

ydQi y dn 



( — ) 



+ ( — ) 
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where (^) et (^) are the cross sections associated 
with incident electrons having helicity plus and minus re- 
spectively. The plus (minus) helicity corresponds to the 
spin of the electron being aligned and in the same direction 
as (opposite to) its momentum. Calculation of the cross 
section requires the knowledge of the amplitudes which 
are derived from the currents in the Feynman formalism. 

The elastic scattering amplitude has two components 
corresponding to the electromagnetic part and to the 
weak part Mz 



M 



(jt',p\he',h^,,k,p, he,h^^ 
= ^Mt(k',p',he',h^,,k,p,he,h^ 



(2) 



where k and k' are the incident and scattered electron, 
p and p' are the incident and recoil proton momentum 
respectively, h^. , and he' , h^, are the electron and proton 
helicity in the initial and final state. 
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where is the Dirac leptonic electromagnetic current: 



u{k' , he') 7" u{k, he) 



(4) 



and J I, em is the hadronic part of the electromagnetic cur- 
rent. The weak amplitude is given by: 



1 
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and the weak currents are obtained from: 



j'^= gl. U{k' ,he') r <k,he) (6) 
J^=g:,u{k\he')rj' U{k,he) (7) 



J 



(8) 



where G is the Fermi constant, gy and g'^ are the weak 
vector and axial charges respectively. For electron scatter- 
ing and at tree level they reduce to = — 1 + 4sin^ 6*,^, 

and g\ — 1 respectively. 

p p 

J I, nc and J n „c5 are the hadronic weak currents. The 
hadronic structure is parametrized in terms of form fac- 
tors: 



F'iQ^) 
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(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



where x — p,n represents a proton p or a neutron n and 
Ux and Ux' are the Dirac spinors for the nucleon in the 
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entrance and exit channel respectively. and are the 

electromagnetic form factors, F^ and F^ are the neutral 

weak vector form factors and and are respectively 
the axial and pseudo-scalar form factors. The latter enters 
in the cross section and asymmetry through aA4z squared 
term which is totally negligible. 

The observables are usually expressed in terms of the 
Sachs form factors , , u and u rather than the 
Fermi and Dirac form factors F , F , F , F : 

1 ' 2 ' 1 ' 2 



GliQ') ^ f;- (Q^) - tF; {Q') 

GI{Q^)^f:{Q^) + f:{Q^) 
G:iQ')^F;iQ^)-rF:m 

g:,{q')^f:{q')+f:{q') 



(15) 
(16) 
(17) 
(18) 



where r is a kinematic factor defined as r = 2 and Mp 
is the proton mass. 

The helicity-correlated cross section is given by: 



< X \UfJ^Uf\x >E 



4'(i^/'^(Q')7M 
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where Uf and uj are the quarks fields. The pseudo-scalar 

form factors G^ being ignored, there are 18 form factors 
to be evaluated: 9 for the proton and 9 for the neutron. 
In order to reduce that number we use charge symmetry, 
assuming that the p and the n are members of a perfect 
isospin doublet. Omitting the Q^-dependence: 



F = F 



F 
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where the summation is performed over all the spin 
variables except the incident electron helicity. The asym- 
metry can then be calculated from Eq. |^: 



GfQ^ 
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in which the dependence has been omitted for clarity 
of notation, e, e' are kinematic factors given by: 



1 



1 + 2(1 + T)tan22|i 
e'- v/^(l + r)(l-e2) 



(22) 
(23) 



9'^ is the electron scattering angle and 6w is the Weinberg 
angle. 

The ultimate purpose of these experiments being to deter- 
mine the strange content of the nuclcon, one must isolate 
the contribution of the s-quark in the nucleon form fac- 
tors. In order to do that, we decompose the electromag- 
netic, neutral and axial currents according to the different 
flavor contributions f = u,d, s: 



F =F 



G\ 



F = F 
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(31) 
(32) 
(33) 



After some algebra [26] , the asymmetry can be finally 
expressed in terms of the electromagnetic, axial and strange 
nucleon form factors: 



GfQ^ 



(l-4sinX.)(l + 



GfQ2 
GfQ^ 



eiGl 



(! + <') 



sg: 



siGlf + riGlf 



rG 



e{Glf + r{G:? 



47raV2 
GfQ^ {l~^siT?e,,)e'Gl 



g: 



G. 



g: 



(34) 



The coefficients i?^ , and are electroweak radiative 
corrections parameters which can be calculated within the 
Standard Model P?] . 

When recoil protons are detected instead of scattered 
electrons, the helicity-correlated cross section becomes: 
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The asymmetry calculation follows the same steps as for 
scattered electron detection. 



3 QED radiative corrections 

3.1 Parity-violating experiment representation in 
leptonic variables 





The aim of the procedure is to get a differential cross 
section d'^a / dQe'dEf.' without any singularity in the full 
electron spectrum. We divide the scattered electron en- 
ergy interval into two regions. The first one, E,,' min < 



Ef>' ^ Ep' , 



where Egi , 



Epi , 



AEgi corresponds 



to "hard photons" with a minimum energy E^' min which 
may be of the order of few MeV. The second one is de- 
fined by Egi cut < Ee' < Egi eias which corresponds to the 
"soft photon" region. The maximum energy of the outgo- 
ing electron corresponds to the elastic peak and is denoted 
by Ee' eias ■ The first requirement is that the integral 



d^a 



df2p'dE„i 



dEp, = 



d^a 



d^a 



dflp'dEpi 



dEp, 



dQpi dEpi 



dEp, 



(36) 



should be as much as possible independent of the cutoff 
energy AEp' . 

The three-dimensional differential cross section which 
appears in the first term in the right hand side of (|36|l is 
obtained from the five-dimensional differential cross sec- 
tion: 



d^<j 



dfipi dEp/ 



dQp,dEp,dQ^ ^^^^ 



corresponding to the bremsstrahlung process e + p — > 
e-|-p + 7. The two Feynman diagrams describing this pro- 
cess are displayed in Fig. ^ The integral defined in j57|l 
may be calculated in the peaking approximation when the 
scattering angle of the detected electron is not too small, 
which is the case for most experiments. In particular, this 
approximation is very good for the PV-A4 experiment at 
forward angle (30° < 6* < 40°) and for the PV-A4 and G° 
experiments at backward angles. The final result is found 
to be jHIEHI: 



Fig. 1. Real-photon emission in virtual-photon exchange 



d^a 



dfipi dEpi 



d^<j 



dQpi dEpi 



peaking 



+ICp 



d^a„{Ep,,Ql) 

dQpi 
d^a^{Ep,Ql^ 



dQpi 



(38) 



where the index stands for the Born elastic differential 
cross section. The term with the s (respectively p) index 
represents the contribution of the left (right) Feynman di- 
agram of Fig. n The coefficients JCg and JCp are kinematic 
factors. 

The second term in the right hand side of H36|l is usu- 
ally expressed as 



d^n 



df2p'dEpi 



dEp, = [l + 5{AEp,) 



d^a„{Ep,Q^) 



dilp 



(39) 



where the theoretical expression oi5{AEei ) is given in j21 
I28| . The numerical calculation of the right hand side of 
Eq. (|3fi|l is shown in Fig. [3 for the PV-A4 parity- violating 
experiment. The minimum value of AEp, in this kinemat- 
ics is about 2 MeV. 

An analytic expression for the three-dimensional dif- 
ferential cross section in the energy range Ep, cut < Ep, < 
Ep' eias can be defined as: 



KdQp/dEpJ 



dQp'dEpJ anal 



a,XOe') + a,{ep,)Ep, + a,{ep,)El 



(40) 
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i{AE,..: 




120 140 160 180 200 

ADC channel 

Fig. 3. Comparison of the experimental (red solid line) and 
simulated spectra (filled grey areas) for the PV-A4 experiment 
|30| . The light-grey area corresponds to the simulated elastic 
ep scattering plus radiative corrections. 



Ec' ela: 



JCs(Ts + /CnC; 



Ee' cut ^ — ^e' elas 



(44) 



^ ^^p^p 

where ai = <Pa^{Ef.^s,Q^)/d^e' , i = s,p. The asymme- 
tries As and Ap are the Born asymmetries calculated for 
the kinematics of the s and p channels through the rela- 
The three parameters a„ , and are fixed using the tions given in the previous section, 
three conditions: 



Fig. 2. Dependence of the right hand side of Eq. I|36|l with 
AE^, in A4 experiment: £^=0.855 GeV, 61^/ =35° 
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(42) 



1 + 5{AEc 



d^a,{Ec,Q^ 
dfipi 



(43) 



3.2 Parity-violating experiment in proton variables 

We describe here the method developed to take into ac- 
count the internal radiative corrections when the proton, 
instead of the electron, is detected. Again, we will ob- 
tain for the proton spectrum a differential cross section 
d'^a / dfipidEpi without any singularity. The extension of 
the method derived for the electrons will give also the 
parity-violating asymmetry in the proton channel. As in 



the electron case, we define Epi cut 
we require the integral l{AEpi) 



E, 



d^<J 



dfipi dEpi 



dEp, 



'p' elas 



d\ 



df2p> dEpi 



AE„' and 



dEp, 



Full simulations performed with the Monte-Carlo method 
and the experimental setup in the angular range between 
40° < 6*° < 30 at Ec = 0.855 GeV pifeO] have shown that 
the final spectrum is, within the experimental resolution, 
independent of the cutoff parameter AEe' when its value 
is increased by a factor 2 to 4. The good agreement be- 
tween the model and the PV-A4 experiment can be seen 
in Fig. 01 

The simulated parity-violating asymmetry is then de- 
fined as: 



d^a 
df2p> dEp' 



■dEr,' 



(45) 



to be as much as possible independent of the energy cut- 
off AEp' . We have to modify the method of the previous 
section for two reasons. First, as we detect the proton, the 
differential cross section is now given by 



dh 



df2p' dEp' 
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d^c 



dQpi dEpi dQj 



dQ-y 



(46) 
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Tll^em = Hk',K') r {^.+ rn,)Y u{k,h,) 



(53) 



while the amplitudes II and IV due to the exchange of the 
Z° with the propagator i{-g^,,, + q^^q^, / M^'^) / {q'^ - M^^) 
have each two different contributions: 




k' k 




P P 




IV 



Fig. 4. Real-photon emission in virtual-Z" exchange. 

Very forward angles of the outgoing electrons are allowed 
when the integration over all the directions of the pho- 
ton is performed, so the cross section has to be calculated 
at the amplitude level to be sure that gauge invariance is 
respected. Secondly, as we are interested to correct the ex- 
perimental asymmetry from the internal radiative contri- 
bution, we need to introduce two more Feynman diagrams 
in the calculation, as shown in Fig. 0] 

The amplitude of the reaction e + p — > e + p -I- 7 is 
written as 



=1, IV 



{k\p\ he'., h^,,pj, h^,k,p,he,\^ (47) 



The four-vectors of the exchanged photon and propa- 
gators are expressed in terms of the kinematic variables: 



q = k — k' — p-y = p' — p 



(48) 



k 



k' + pj ^ k + p - p' (49) 



The amplitudes I and III due to the exchanged photon 
with the propagator —ig^u' j q^ have one term: 
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In the energy range where the parity-violating exper- 
iments are performed (0.1 < < l{GeV/c)^ , terms 
proportional to l/M^^ are neglected, therefore: 



Mil ~ ie 
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and 
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(63) 



The total amplitude is the sum of two terms 
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The interference of these two terms will produce the parity- 
violating asymmetry. The Parity Conserving amplitude 

PC 

M. is due to photon exchange and it contains a part of 

PV 

the exchange. The Parity Violating amplitude M is 
due to part of the Z° exchange contribution in the Feyn- 
man diagrams II and IV. Explicitly, this amplitudes is 



M 



PV 
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PV 
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Now we are able to calculate the integral 



^p' cut 



^^dE, 

df2j)' dEn' 



d-^a 



df2p' dEp' 
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PV 
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M 
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IV 
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G _ 
2^/2 x\ 



[J.ncTll\ + JuncTi;^,,) e;{p^,h,) (67) 

The differential cross section is then calculated in the lab- 
oratory system in terms of the amplitudes by: 



dfipi dEpi df2~f 



32(27r) 

1 

4 



M\k\ 



\E^ 



Eg' + Ury ■ — k) 



(68) 



where the summation is performed over all the helicity 
states of the incoming electron, the target, the outgoing 
proton and the outgoing photon. The differential cross sec- 
tion of the outgoing proton is then expressed as in H46(l . 
after integration over all the photon angles. 

The parity-violating asymmetry is calculated in a sim- 
ilar way. First we calculate the differential cross section as 
a function of the beam helicity hg = ±1/2 : 



df2p' dEp' df^'j 



32(27r)5 



1 



2 \E^ 



\P'\E, 
M\k\ 



- • (p' — fc) 



(69) 



The prime index over the summation means that the sum 
is performed over all the spin variables except the incident 
electron helicity. The parity-violating asymmetry of the 
proton spectrum then reads: 



for any value of AEp' 0. As in the electron case, the 
integral in the energy range Ep^ cut < Ep^ < Ep> eias is 
proportional to the Born elastic differential cross section: 



''^'^^ dE.^AiAEp,)''^^^^-'^'^ 



dfin' dErj' 

p' cut ^ ^ 



dQn 



(73) 



Its calculation is given by the following ratio 



A{AEp>) 



ith inn 



K{AEp,) 

dQ^ 

dQpi dEpi dfi-y 



dn^ 

df2p' dEp' df2'^ 



K-iAE^^ 
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'P'j 
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/2f 



(74) 



(75) 



The meaning of K.{AEpi) is clear. For each value of 9p 



Epi , 9.y and i 



, the value of AEp. 



is equal to E, 



p' elas ' 



Epi. The three body kinematics give the energy of the 
photon E^ and the complete kinematics of the outgoing 
electron 9^' , 0e' and E,,/ through the energy-momentum 
conservation. Comparison with the elastic scattering e -|- 
p — > e -t- p reaction at the same angle gives the value of 
AEg, = Eg, das - Eg,. The ratio K. = e*"-'-+*«/(l - 
S vacuum/ '2')'^ is the attcnuation factor, which depends on 
AEg, , induced by the internal radiative correction on the 
electron side. It is a generalization to all orders of the 
1 + S term of Eq. (|43|l as it can be seen if we make a Taylor 
expansion of /C. Finally the attenuation factor A{AEp, ) as 
defined in equation (|74|l is the average attenuation factor 
when we integrate over all the directions of the photon. 
The explicit formulae used in the code to calculate Sygrtex , 
Sn and Syacuum are taken from CTT] : 



Or = - 

TT 



/(AE^ 
V EgEgi 



\mlJ J 2 \Eg,J 



8 



J. Arvieux et al.: Electromagnetic radiative corrections in parity- violating electron-proton scattering 



X 10 — 
0.672 - 



0.671 




0.668 



0.666 



0.314 0.315 0.316 0.317 0.3U 

Tp (GeV) 



Fig. 5. T{AEpi) as a function of the kinetic energy of the 
detected proton for 6pi =60° 
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(76) 

(77) 
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(79) 



The value of the integral I{AEpi ) as a function of AEpi 
has been performed for 48° < 9p' < 77°. It is plotted in 
Fig- El for one scattering angle of the detected proton. The 
value of the cutoff parameter is chosen so that this integral 
reaches its minimum value. 

As in the electron case, we assume that for the kinetic 
energy range of the scattered proton Tpi eias — AEpi < 

p' elas 



Tpi < T, 



^ df2p' dEpi / anal 



^o(^P') + '^li^p') {Tpi — Tp' elas) 



+ a^{9pi) {Tpi —Tpielas)^ 



(80) 



The determination of the three coefficients a^, a-^ and 
is obtained by the following conditions: 



i) at Epi — Epi cut ■ 
\dQ„idEwJ V 



df2p' dEp/ ) anal 



(81) 




Tp (GeV) 

Fig. 6. Parity-violating asymmetry as a function of the ki- 
netic energy of the detected proton for 6pi=60°. The star * 
represents the Born asymmetry. 



at Epi — Epi cut ' 
d / d^a 



KdQwdE^i) ^ dE„i V 



dEpi V df2pi dEpi 



7r , \df2r>'dE„i J anal 

p cut ' " 



dEp! \df2p' dEp' / anal 



(82) 



III] 



d^a„iEe,Q^) 
df2„i 



(83) 



The parity-violating asymmetry is calculated through 
the relation H7()() for Epi < Epi cut and near the elastic 
peak, its value is linearly interpolated between its value at 
Epi — Epi cut and the Born asymmetry calculated at Epi = 
Epi elas ■ The variation of this asymmetry as a function of 
the kinetic energy of the scattered proton is plotted in 
Fig. Elfor one angle. 



4 Application to the experiment at 
forward angles 

The experiment performed in Hall C at Jefferson 
Lab, measures the parity-violating elastic electron scat- 
tering from the nucleon. Asymmetries of the order of one 
part per million from scattering of a polarized electron 
beam are determined using a dedicated apparatus. It con- 
sists of specialized beam monitoring and control systems, a 
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cryogenic hydrogen target and a superconducting, toroidal 
magnetic spectrometer equipped vifith plastic scintillation 
counters as well as fast readout electronics for the mea- 
surement of individual events. 

In the forward-angle configuration, a polarized elec- 
tron beam of 40 fiA with an energy of 3.031 ± 0.001 GeV 
was used over the measurement period of 700 h. It was 
generated by a strained GaAs polarized source with 
32-ns pulse timing (rather than the standard 2 ns) to allow 
for time-of-flight (tof) measurements. The average beam 
polarization, measured with a M0ller polarimeter j34j in 
interleaved runs, was 73.7 ± 1.0%. Helicity-correlated cur- 
rent and position changes were corrected with active feed- 
back to levels of about 0.3 parts-per-million (ppm) and 10 
nm, respectively. Corrections to the measured asymme- 
try were applied via linear regression for residual helicity- 
correlated beam current, position, angle and energy vari- 
ations, and amounted to a negligible total of 0.02 ppm; 
the largest correction was 0.01 ppm for helicity-correlated 
current variation. 

The polarized electrons scattered from a 20 cm liq- 
uid hydrogen target; the recoiling elastic protons were 
detected to allow simultaneous measurement of the wide 
range of momentum transfer 0.12 < < 1.0{GeV/c)'^ . 
This was achieved using a novel toroidal spectrometer de- 
signed to measure the entire range with a single field set- 
ting and with precision comparable to previous experi- 
ments. The spectrometer included an eight-coil supercon- 
ducting magnet and eight sets (or octants) of scintillat- 
ing detectors. Four octants (numbered 1-3-5-7) and their 
associated electronics were built by the North- American 
(USA, Canada) part of the G° collaboration and four oc- 
tants (2-4-6-8) and their associated electronics were built 
by the French (IPN Orsay, LPSC Grenoble) part of the 
G° collaboration. Each set consisted of 16 scintillator pairs 
used in coincidence to cover the range of momentum trans- 
fers (smallest detector number corresponding to the lowest 
momentum transfer). The scattering angle varies from 52 
to 76 degrees, depending on detector number. Because of 
the correlation between the momentum and scattering an- 
gle of the elastic protons (higher momentum corresponds 
to more forward proton scattering angles), detector 15 
covered the range of momentum transfers between 0.44 
and 0.88 (GeV/c)^ , which we divided into three tof bins 
with average momentum transfers of 0.51, 0.63 and 0.79 
(GeV/c)^ . For the same reason, detector 14 had two elas- 
tic peaks separated in tof with momentum transfers of 0.41 
and 1.0 {GeV/c)'^ ; detector 16, used to determine back- 
grounds, had no elastic acceptance. Custom time-encoding 
electronics sorted detector events by tof; elastic protons 
arrived about 20 ns after the passage of the electron bunch 
through the target. A typical time-of-flight spectrum is 
shown in Fig. [7| The spectrometer field integral and ulti- 
mately the calibration {AQ'^/Q^ = 1%) was fine-tuned 
using the measured tof diff'erence between pions and elas- 
tic protons for each detector. All rates were corrected for 
dead-times of 10 — 15% on the basis of the measured yield 
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Fig. 7. Experimental yield as a function of tof for detector 8. 
The elastic protons correspond to the rightmost peak. 

dependence on beam current; the corresponding uncer- 
tainty in the asymmetry is 0.05 ppm. The final results 
of the G'^ forward-angle experiment are shown in !1 II . 

Radiative corrections for G° have been estimated in a 
simulation using the GO-GEANT package j^H- The elec- 
tron can, in principle, loose all its energy through radia- 
tion, but the probability that it looses 500 MeV or less is 
96%. Moreover, 60% loose 1 keV or less. The few events 
for which the electron energy loss is more than 500 MeV 
correspond to proton having times of fiight out of the G° 
experimental cuts, thus they are not considered in our 
calculation. Therefore the GEANT simulations have been 
done in the energy interval Ei„c= 2.5-3.0 GeV only and 
for recoil proton angle 9p' = 48° — 77° and energy Tp/ =2 
MeV-Tp5 . The cross sections have been interpolated for 
intermediate values using a spline method. In order to ob- 
tain rates, each event (number j) is normalized through a 
weight Wj proportional to the cross section 



[Emax{Ej"''^ , 9j) — EjniniEj"''' , 9j)] (84) 

where £ is the luminosity, A(l> is the polar angle opening 
and Nt is the number of drawings. In the case of clastic 
scattering (Born term), the weight is simply given by 



5 Results 

5.1 Time-of-flight spectra 

Two calculations are performed without and with RC: 

- in the first case, the time-of-flight of elastically 
scattered protons, without any energy loss nor radiative 
corrections is calculated. The width of the peak is essen- 
tially given by the experimental resolution as calculated 
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Fig. 8. Simulated elastic-proton time-of-flight distributions 
for Det. 8 (top: FR, bottom: NA). The gaussian fits are per- 
formed to extract the position of the mean. 



Fig. 9. Comparison of time-of-flight spectra (FR Det. 8) with- 
out (black dash-dotted line) and with (red/grey solid line) ra- 
diative corrections. 



by GEANT. The elastic peak is fitted with a Gaussian al- 
lowing to determine the position of the maximum, in order 
to define cuts within which the asymmetry is calculated. 
The resulting spectra are shown in Fig. (HJ where only de- 
tector 8, corresponding to the middle of the focal plane, is 
shown for reference. The only difference in the two spec- 
tra is a binning of 250 ps for the French (FR) electronics 
(top) and a binning of 1 ns for the North- American (NA) 
electronics (bottom). 

- in the second case, the proton tof spectra are 
calculated after applying energy losses and full RC. The 
result obtained for detector 8 is shown in red (grey) in 
Fig. El overlaid to the pure elastic spectra (in black). 
One should notice the following paradox: Inelastic protons 
have, by definition, an energy lower than elastic protons 
and therefore, a smaller velocity. They should then appear 
on the right side of the elastic proton peak, correspond- 
ing to a longer time-of-flight. In fact, due to the effect of 
the magnetic field and to the geometry of the collima- 
tors, inelastic protons have a shorter trajectory and they 
reach a given detector faster than the elastic ones. This is 
confirmed in the experimental data. 



5.2 Asymmetries 

The asymmetry is calculated using Eq. (|34|l . with the fol- 
lowing numerical values 0: 
sin^e,,, = 0.23117 



Gf = 1.16639 X 10-5 

Rl -0.0520, Rl = -0.0143, R^"^ ^ 0. 

and G^^ are parametrized with a dipole form accord- 
ing to |2| ■ A discussion of the latest electromagnetic form 
factors can be found in jjl]. The strangeness content pa- 
rameters are from Hammer et al. |36) with /is = —0.24 
and Ps = —2.93. These values have been taken from a re- 
view paper by Kumar and Souder 0. These parameters 
are used here only as an example of a strange asymme- 
try calculation. Electromagnetic radiative corrections are 
rather insensitive to the clectroweak parameters. Detailed 
calculations will be shown on the FR detectors spectra 
only. Fig. EH shows the asymmetry distribution without 
(black line) and with RC (red/grey line). In the elastic 
case, the asymmetry reduces to the one calculated from 
the Born term only, and it can be compared directly to 
theoretical models. 

The mean asymmetry value is plotted in Fig.^Jfor the FR 
detectors. The effect of radiative corrections is to increase 
the average asymmetry, following the increase in . The 
ratios between elastic and RC corrected asymmetries are 
given in the following Tabled 

The 'ideal' procedure to analyze the data would be to 
calculate an experimental asymmetry from the data after 
removing all background, leaving an experimental peak 
including radiative emission, and then to multiply the cor- 
responding asymmetry by the ratio Apj j A^c given below. 
One problem is that, when removing background using a 
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Fig. 11. Mean value of the asymmetry for each detector. The 
black stars represent the elastic case, the red (grey) triangles 
the radiative one. The "french cuts" on the tof distributions 
have been applied. 



Fig. 10. Asymmetry distributions (in ppm), for each detector. 
The black dash-dotted line represent the elastic case, the red 
(grey) solid one the radiative case. 

pure fitting procedure, one removes also part of the inelas- 
tic tail and therefore, by using the above procedure, one 
will overestimate the radiative corrections. A quantitative 
estimate of this effect is given below. 
The asymmetry increase is of the order of 0.5-1.0 % for 
detectors 1-9, reaching 2.0 % for detector 12% and up to 
3.0% for detector 14. These ratios should be almost inde- 
pendent of the model chosen and therefore valid for the 
no-strangeness value Aq. It is not clear if the dispersion 
between correction factors between 2 adjacent detectors 
(e.g. between Det. 8-9-10 or 10-11-12), which is of the 
order of 0.3 %, is an indication of the present statisti- 
cal/systematical errors or if it is a genuine effect due to 
differences in acceptance. 



5.3 Uncertainty estimate 

An error estimate is made based on the assumption that 
the elastic cuts have a 10% uncertainty. Therefore the ra- 
diative corrections are calculated for cuts which are 5% 



larger than the elastic cuts (by increasing the upper limit 
by 2.5% and decreasing the lower limit by 2.5%) and 5% 
smaller than the elastic cuts (by decreasing the upper limit 
by 2.5% and increasing the lower limit by 2.5%). Then we 
take the ratio of these two quantities for each detector. 
This should represent an upper limit of the radiative cor- 
rection uncertainties since the elastic cuts are known to 
better than 10%. The corresponding uncertainty would 
vary slowly from 0.1% for Det. 1 to 0.5% for Det. 13, 1% 
for Det. 14 and between 0.0% and 0.7% for Det. 15, de- 
pending on the cut. An alternative error estimate is 
obtained by making a global fit of the ratio Aei/Aj^c with 
a polynomial and assuming that the difference with the 
actual RC correction is due to systematics: in that case 
the uncertainty is globally estimated to be of the order 
of 0.1-0.3% or 10% of the actual correction depending on 
detector number. 

Another problem which has been investigated is the 
one of correction double-counting. If the background un- 
der the elastic peak is removed by a pure fitting procedure, 
it will also contain the RC tail contribution to the peak. 
Therefore the corresponding elastic asymmetry should not 
be corrected for RC effects. In order to estimate the sen- 
sitivity of the RC corrections, at the border of the elastic 
peak, we have calculated the RC corrections by adding or 
removing 1 ns from the elastic cuts. This effect has been 
estimated to be about 2% of the RC corrections which are 
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Table 1. Ratio of asymmetries Aei/Anc as a function of de- 
tector number where A^i is the elastic (Born term) asymmetry 
and Arc is the asymmetry corrected for radiative emission. 
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themselves of the order of 2%, so that double counting can 
be neglected at first order. 



6 Summary and conclusions 

We have calculated the full electromagnetic radiative cor- 
rections for elastic ep scattering in leptonic or hadronic 
variables; a performable code has been constructed to ex- 
tract the parity-violating asymmetry from the experimen- 
tal measured asymmetry. The comparison between the 
simulation results and the data in the kinematic configu- 
ration of the PV-A4 experiment validates our procedure. 
Radiative corrections for the GO parity-violating elastic 
scattering experiment have been estimated by feeding our 
model calculations through a Monte Carlo detector sim- 
ulation. This code could also be used for the next asym- 
metry measurement in the backward-angle configuration 
of GO. 
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